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ABSTRACT 

Observations of four comets made with the Far Ultraviolet Spectroscopic Ex- 
plorer show the rotational envelope of the (0,0) band of the CO Hopfield-Birge 
system (C l Y> + - A 1 S+) at 1088 A to consist of both "cold" and "hot" com- 
ponents, the "cold" component accounting for ~75% of the flux and with a 
rotational temperature in the range 55-75 K. We identify the "hot" component 
as coming from the dissociation of CO2 into rotationally "hot" CO, with electron 
impact dissociation probably dominant over photodissociation near the nucleus. 
An additional weak, broad satellite band is seen centered near the position of 
the P(40) line that we attribute to CO fluorescence from a non-thermal high 
J rotational population produced by photodissociation of formaldehyde into CO 
and H2. This process also leaves the H2 preferentially populated in excited vibra- 
tional levels which are identified by fluorescent H2 lines in the spectrum excited 
by solar O VI A1031.9 and solar Lyman-a. The amount of H 2 produced by H 2 CO 
dissociation is comparable to the amount produced by photodissociation of H 2 0. 
Electron impact excitation of CO, rather than resonance fluorescence, appears 
to be the primary source of the observed S 1 E + - X 1 E + (0,0) band at 1151 A. 

Subject headings: comets: individual (C1999 Tl, C/2001 A2, C/2000 WM1, 
C/2001 Q4) — ultraviolet: solar system 
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INTRODUCTION 



In previous papers ( Feldman et al.ll2002l ; IWeaver et al.ll2002uFeldmanll2005l ). we reported 
on the spectra of four comets observed by the Far Ultraviolet Spectroscopic Explorer (FUSE). 
Launched in June 1999, FUSE operated through October 2007 providing an orbiting capa- 
bility with a spectral resolution better than 0.4 A in the wavelength range 905-1187 A 
together with very high sensitivity to weak emissions making possible both the search for 
minor coma species and extensive temperature and density diagnostics for the dominant 
species. Our initial reports on the three comets observed in 2001 focused on three Hopfield- 
Birge band systems of CO not previously observed in comets, the B l Y, + - X 1 S + (0,0) band 
at 1150.5 A, the C X T, + - X 1 S+(0,0) band at 1087.9 A, and the weak E 1 !! - X 1 S+(0,0) 
band at 1076.1 A, three lines of H2 fluorescently pumped by solar Lyman-/? radiation, and 
upper limits to emission from O VI and Ar I. The rotational envelopes of the CO bands 
are resolved and appear to consist of both cold and warm components, the cold component, 
with a rotational temperature in the range 55-75 K, accounting for over 75% of the flux 
and presumably due to native CO released by the nucleus. The warm (~500 K) component 
suggested a CO2 photodissociation source but this was considered unlikely because the low 
rate for this process in the solar radiation field would have required a very high relative 
abundance of C0 2 . 

O VI emission at 1031.9 and 1037.6 A was searched for but only the 1031.9 A line 
was marginally detected in C/2000 WM1 at about the level predicted by a model used to 
explain the X-ray emission from co mets as being produced by char ge exchange between solar 
wind ions and cometary neutrals (IKharchenko fc Dalgarnoll200lf ). The higher signal/noise 
ratio in the spectrum of comet C/2001 Q4, observed in April 2004, enabled an accurate 
determination of the wavelength of this feature which suggested that it was the H 2 (1,1) Q(3) 
Werner (C 1 n n — X line, which can be fluorescently pumped by either solar N III A989.8 
or O VI A1031.9 . The latter requires the H2 to be in the first vibrationally excited state 
Liu fc Dakarnolll996h . The Q(3) lines of the (1,3) band at 1119.1 A and the (1,4) band 
at 1163.8 A are also detected as would be expected for fluorescent excitation. However, 
the ground state H2 abundance in the coma, estimated from the solar Lyman-/? pumped 
fluorescence of the Lyman (B X E+ — X X S+) (6,1) PI line, is insufficient for excitation by 
solar N III (a factor of 10 weaker than O VI), which implies the formation of H2 in the v = 1 
state and excitation by solar O VI. 

The molecular channel of the photodissociation of formaldehyde (i.e., the channel lead- 
ing to H 2 and CO products) has been well studied both experimentally and theoretically. Its 
products are characterized by a non-thermal rota tionally "hot" CO molecule and a vibra- 
tionally excited H 2 molecule (Ivan Zee et al.lll993l ). The CO rotational distribution can be 
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represented by a gaussian centered near J = 40 with a half-width at half maximum in J of 
about 20. Fluorescence from such levels in the P-branch (the R-branch is shortward of the 
detector wavelength cut-off) is clearly seen in two of the comets observed by FUSE and we 
propose that this feature, together with the vibrationally excited H 2 , constitute an unam- 
biguous signature of H2CO in a cometary coma. Formaldeh yde has been observed in these 



comets in both the inf rared and millimeter spectral regions (IBiver et al.l 120061 ; iMilam et al. 



20061 ; iGibb et al.ll2007l ). The capability of FUSE to identify H 2 CO in a cometary coma pro- 
vides a new means to determine its production rate relative to that of water and to compare 
these results with ground-based observations in the infrared and millimeter regions of the 
spectrum. 

Since the thermal "hot" component of the observed CO emission is not produced by 
H 2 CO dissociation, we return to C0 2 as a possible parent, particularly as the rotational 
temperatur e is what would be expected based solely on angular momentum conservation 
arguments (IMumma et al.lll975l ). We suggest that within the FUSE field-of-view, photo- 
electron impact dissociation of C0 2 has a comparable rate to photodissociation and we base 
this conclusion on estimates of the mean photoelectron flux derived from the other prominent 
CO band observed in these spectra, the B — X (0,0) band at 1151 A. We thus have unique 
diagnostics in which the shapes of these molecular bands can be used to infer the column 
abundance of CO from both its native source as well as from two of its direct molecular 
parents, together with the photoelectron contributions to excitation and dissociation in the 
inner coma. 



2. OBSERVATIONS 

Three long period comets were observed during 2001, C/1999 Tl (McNaught-Hartley), 
C/2001 A2 (LINEAR), and C/2000 WM1 (LINEAR), with C/2001 A2 being the most active 
of the three and showing the brightest far-ultraviolet emissions. In December 2001, two of 
the reaction wheels on the spacecraft failed resulting in a shutdown of observatory operations 
until a new pointing control program could be devised and implemented a few months later. 
However, no moving target observations were attempted until early in 2004 when comet 
C/2001 Q4 (NEAT) was observed. The observation parameters for all four comets are 
summarized in Table [TJ 

Observations of each comet were made over contiguous spacecraft orbits during which 
the center of the 30" x 30" spectrograph aperture was placed at the ephemeris position of 
the comet and tracked at the predicted sky rate. Due to the extended, though not uniform, 
nature of the cometary emissions within the aperture, the effective spectral resolution is 
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~0.25 A. For each comet, all of the orbits of data were co-added and the extracted fluxes 
were converted to average brightness (in rayleighs) in the 30" x 30" aperture. In addition, the 
data were separately co-added for only the time during which the spacecraft was in Earth 
shadow. These "night-only" spectra make it possible to differentiate cometary emissions of 
H I, O I, N I, and, in se cond order, He I, fro m the same emissions produced in the daytime 
terrestrial atmosphere (IFeldman et al.ll200ll ). For the CO and H2 emissions this is not a 
problem and the entire data set for each comet was used. 

All of the data were rep rocessed using the final version of the FUSE pipeline, Cal- 
FUSE 3.2 (IDixon et al.l 120071 ). Among the enhancements incorporated in this version were 
temperature-dependent wavelength corrections, time-dependent corrections to the detector 
effective areas, and a more complete evaluation of detector background. The latter was par- 
ticularly significant in allowing reliable extractions of the fluxes of weak emission features. 
The overall absolute flux calibration, based on standard stars and adjusted for the date 
of observation, is better than 5%. Table [1] also lists the extracted brightnesses of the CO 
emission features discussed below. 



3. SPECTRA 



3.1. Data 



The CO C — X (0,0) band, from the total exposure, is shown for each comet in Figured! 
The P and R branches appear well separated, characteristic of a rotational temperature of 
55-75 K. There also appear broad wings to the band indicative of a rotational temperature 
in the range of 500 to 600 K. The figures also show a model fit, assuming s olar resonance 
fluorescence, obtained using the procedure described by IFeldman et al.l (120021 ). The derived 
column densities and rotational temperatures are given in Table [21 In the two most active 
comets, C/2001 A2 and C/2001 Q4, there is also an additional broad feature centered at 
~1089.4 A, that we identify as P-branch lines from hi gh-lying J states . The position of the 
P(40) line, calculated using the rotational con stants of lEidelsberg et al.l (ll99ll ). is also shown 
in Fig. [TJ We note that the lEidelsberg et al.l wavelengths have been verified experimentally 
only up to J = 25. The corresponding high J R-lines are off the short wavelength edge of the 
LiF2a detector. A "hot" thermal rotational distribution arises in molecular photodissociation 
in order to conserve angular momentum when the product molecul e has a much smalle r 
moment-of- inertia than the parent molecule as is the case for CO2 (IMumma et al.lll975l ). 
The non-thermal component is identified with CO as a dissociation product of H 2 CO, and 
the column abundance of this component, also given in Table [21 is derived assuming an 
equal contribution from the unobserved R-branch. In this case, an optically thin band 
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fluorescence efficiency is used as the CO population is spread over many rotational levels 
and the absorption is not likely to be optically thick. 



3.2. Photodissociation of Formaldehyde 



As noted above, the photodissociation of formaldehyde has been intensively studied, 
both experimentally and theoretically. The lowest dissociation threshold (~3600 A) leads to 
H 2 + CO (the "molecular" channel) wh ile ~2700 cm" 1 high er (-3300 A), H + HCO (the 



"radical" channel) can also be produced (Ivan Zee et al.lll993l ). The nature of the H2 and CO 
products differs markedly above and below this second energy th resho ld. Consider fi r st the 
region below the radical channel threshold. iBamford et al.l (119851 ) and iDebarre et al.l (119851 ) 
have used laser-induced fluorescence spectroscopy and coherent anti-S tokes Raman sc attering 
to experimentally study the product state populations of CO and H 2 . IBamford et al.l found a 
nonthermal highly excited CO rotational population with a Gaussian distribution centered at 
J = 40 and a half - width at half maxim um of 20 J units, but with little vibrational excitation. 
Later work (e.g.. Ivan Zee et al.lll993l ) showed that the shape was dependent on the energy 



of the absorbed photon. It sh ould be n oted that the a bsorption of H 2 CO in this region is in 
discrete bands (lRogerslll99Cll ). For H2, IDebarre et al.l found that the primary products are 
the ortho ( J = 3 and 5) states of v — 1. For ph otolysis at 3390 A they found the fractional 
population of the v = 1, J = 3 level to be 11%. iButenhoff et al.l (jl990l ) extended this work 
to additional wavelengths and found a peak H 2 population in v = 1, J = 3. 

Above the radical channel threshold, H 2 and CO are found to have markedly different 
populations with the non-translational energy going principally to vib rational excitation o f 
H 2 while the CO r otational state population shifts to lower J states (Ivan Zee et al.lll993l ). 
Zhang et al.l (120051 ) demonstrated that the population of vibrationally excited H 2 , up to 
v = 8 or 9, increased as the energy of the exciting photon increased with a corresponding 
increase in the low-J state population of CO such that the non-translational energy was 
conserved. H owever, the peak a t J = 40 remained and the low-J population was still highly 



non-thermal. IChambreau et al.l (120061 ) showed that the excited vibrational levels of H 2 also 
had peak rotational populations at J = 5 or 7, depending on the photolysis energy. 

In the solar radiation field in the coma of a comet, the H 2 and CO state populations 
derived from H 2 CO dissociation will re flect a superposition o f populations, weighted by the 
solar flux and absorption cross section (IHuebner et al.lll992l ). Such a calculation is beyond 
the scope of the present paper. The broad, weak feature seen longward of the C — X 
(0,0) band in comets C/2001 A2 and C/2001 Q4 in Fig. [1] is centered near the position 
of the P(40) line of this band, and the non-thermal high- J population of this feature is 
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an unambiguous fingerprint of CO produced from H2CO photodissociation. Although the 
spectrum of C/2000 WM1 is significantly noisier, the data also suggest enhanced emission on 
the red wing of the band in this comet. Further evidence for this sou rce is the simul taneous 
detection of fluorescence from vibrationally excited H 2 in this comet (ILiu et al.ll2007l ) as well 
as in comets C/2001 A2 and C/2001 Q4. 



3.3. Fluorescence of Vibrationally Excited H 2 



In Figure [2] the region of the O VI doublet in comet C/2001 Q4 is shown, and the 
H 2 (1,1) Q(3) line is identified, separated by 0.06 A from the position of the O VI line 
at 1031.92 A. The uncertainty in the FUSE wavelength scale is ~0.01 A. The same figure 
also shows t he corresponding (1,3 ) Q(3) and (1,4) Q(3) lines at 1119.09 and 1163.79 A, 
respectively (ILiu fc Dalgarnolll996l ). Note that there is no emission at the expected position 
of the O VI line at 1037.62 A. These lines a re qu ite weak and with the exception of a 3-a 
detection in comet C/2000 WM1 (ILiu et al.l 120071 ) do not appear in the spectra of the two 
other comets observed by FUSE. 

The absence of these lines in the other two comets is the result of the heliocentric 



Doppler shift of the exciting solar line, the Swings effect (IFeldman et al.ll2004l ). rather than 
the lack of H 2 in the coma . This can be illustrated with the solar O VI A1031.9 line profile. 
Doschek fc Feldmanl (120041 ) have used the recent solar spectra taken by the SUMER instru- 
ment on SOHO with an instrumental resolution of ~0.1 A, to determine the width of this 
line under quiet sun conditions to be 0.202 A. This result is from an observation of a small 
part of the solar disk near the center, while for the calculation of cometary excitation rates 
the whole disk flux is more us eful. Fortuitously, w e are able to take advantage of a minor 
operational problem of FUSE (ISahnow et al.ll2000l ). a sensitivity at certain orientations on 
the sky of the silicon carbide (SiC) channels to scattered solar radiation, to determine the 
whole disk line shape. The observation of comet C/2001 Q4 resulted in such an orienta- 
tion and the SiClb and SiCla channels (905-992 and 1005-1090 A, respectively) measured 
a large scattered solar flux in all three apertures. Both the high (HIRS, l."5 x 20") and 
medium resolution (MDRS, 4" x 20") apertures obtained good S/N spectra of the O VI lines 
at instrumental resolutions of 0.035 and 0.075 A, respectively. From both of these channels 
we derived a width of 0.200 A for the solar O VI A1031.9 line, in excellent agreement with 
the result of iDoschek fc Feldmanl . The MDRS profile is shown in Figure El together with 
the Doppler shifted positions of the H 2 (1,1) Q(3) line for each comet. Note that the profile 
shown in the figure suggests that the solar line is not adequately represented by a single 
Gaussian. Clearly, the detection of the (l,v") lines in comet C/2001 Q4 was favored by its 
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negative heliocentric velocity and by the relatively high activity level of this comet at the 
time of observation. 



Wolven et al.l ( 19971 ) noted that solar Lyman-a (1216.57 A) could pump vibrationally 
excited H 2 in the v = 2 level via the Lyman band (1,2) P(5) and R(6) lines, and used this 
mechanism to identify fluorescence from the upper level in Hubble Space Telescope spectra of 
Jupiter following the impact of comet Shoemaker-Levy 9 in July 1994. The same mechanism 
has recently been found to operate in planetary nebulae where the FUSE spectra sh ow the 
(1,1) P(5) and R(6) lines at 1161.82 and 1161.95 A, respectively rtLupu et al.ll2006h . One 
of the previ ously unid e ntifie d lines reported in FUSE spectra of comets at 1161.83 A was 
identified by lLiu et al.l (120071 ) as a blend of these two H2 lines (see Fig. [2]). The vibrationally 
excited H9 would m ost likely come from H 2 CO rather than H 2 photodissocation as pro- 
posed by iLiu et al.l This feature has roughly the same brightness in comets C/2001 A2 and 
C/2001 Q4, and we note that the brightnesses of the CO P(40) featu re in these two comets 
are also comparable. From the branching ratios given by iLupu et al.l . we would also expect 
to see the (1,1) R(3) line at 1148.70 A and the (1,1) P(8) line at 1183.31 A from the same 
excited J levels, and indeed we do. The latter is considerably weaker than the former due 
to the lower population of the even-J (para) states as expected. 



Liu et al.l also identified several other cometary lines as being pumped by solar Lyman-a 



from higher vibrat ional and rotatio nal levels of H2. Almost all of the strong unidentified fea- 
tu res tabula ted by iFeldmanl (120051 ) are amongst them. Our analysis supports the conclusion 
of lLiu et aLl that the vast majority of emissio n features in the FUSE spectra of comets can be 
attributed to H 2 . However, unlike ILiu et al.l . who proposed that the excited vibrational and 
rotational levels of H 2 arise from the photodissociation of H 2 0, we suggest instead that these 
excited H 2 states are produced during the photodissociation of H 2 CO because their presence 
in the cometary coma is consistent wi th the level pop ulati ons measured dur i ng Ho CO 
photodissociation in the laboratory by IZhang et al.l (120051 ) and IChambreau et al.l (120061 ) . 



3.4. Electron Impact Excitation of CO and Dissociation of C0 2 

The CO B — X (0,0) band at 1151 A is weaker than the C — X (0,0) band by a factor of 
3 to 4 and thus recorded with poorer signal / noise in all four comets as shown in Figure HI It 
is surprising that this band is observed at all since its fluorescence efficiency is 44 times lower 
than that of the C — X (0,0) band (the predicted spectra using the CO column densities given 
in Table [2] are also shown in the figure). There must be an additional source of excitation, 
and since the band shape is consistent with the low temperature CO component, that can 
only be electron impact excitation of CO. Photoelectrons in the inner comae of comets have 
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been considered by several authors (e.g., lAshiharalll978 



detected in the Vega fly-by of lP/Halley (IGringauz et al 



^orosmezey et al.lll987l ). and were 



19861 ). but their role in excitation, 



dissociation and ionization is difficult to quantify 

An equivalent to the fluorescence efficiency, or g-f actor, may be defined by an electron 
excitation rate per molecule per second, g'\ 



g ' = J J F{E)a{E)dEdVt 



where F(E) is the photoelectron flux in electrons s~ 1 cm~ 2 sr~ 1 eV _1 and a(E) is the cross 
section in cm 2 . In general, the photoelectron flux is a function of cometocentric distance and 
is anisotropic so that g' will vary over the spectrograph field-of-view. 

Rather than attempt to use modeled photoelectron fluxes, we instead use the observed 
CO B — X (0,0) band brightnesses to derive an effective g' over the FUSE aperture but 
then use the modeled electron energy distribution to determine the relative excitation rates 
for other processes. In particular, we are interested in the electron impact dissociation rate 
of CO2. Relevant cross sections in analytic form fo r electron impact on CO and CO2 are 
available from the compilation of IShirai et al.l (120011 ) . Figure [5] shows the cross sections for 
electron impact excitation on CO of the two bands under discussion and that for electron im- 
pact dissociation of CO? into the met a st able a 3 II state of CO. A typical photoelectron energy 
dependence, from iKorosmezey et al.l (119871 ). is also shown. The a 3 IT state is t he upper level 



of the CO Cameron bands, whose excitation in comets has been dis cussed by IWeaver et al. 



(1l994f ) . The total CO2 dissociation cross section has been studied by I Cosby fc Helml (119921 ). 
These authors note that the branching ratio into the a 3 IT state o f CO is 27% of the total 
electron dissociation rate, the same as for solar photodissociation (IHuebner et al.lll992l ). 



The derived values of g', reduced to 1 AU, are given in Table [3j For comparison, the 
optically thin solar fluorescence efficiencies are also given. Figure [5] shows that the steep 
increase of photoelectron flux towards lower energies favors the B — X transitions over the 
C — X. Even so, the electron impact contribution to C state excitation will reduce the 
derived CO column abundances by 5-15% and this is reflected in the values given in Table [2j 
These numbers remain uncertain as both the photoelectron energy distribution and flux will 
vary with increasing column density towards the nucleus. This will also affect the ratio of 
total C0 2 dissociation to that fraction leading to the a 3 II state of CO. The g' values for these 
processes, for the photoelectron model adopted, are also given in Ta ble El For comparison 
with photodissociation at 1 AU, the solar maximum value given by IHuebner et al.l (119921 ) 



is 2.56 x 10 6 molecule 1 s 1 , so electron impact dissociation and photodissociation are 



comparable. 
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Evidence for a significant role of photoelectrons in the inner coma has often been noted 
in the literature. Two particular cases in which photodissociation and photoionization were 
inadequate to account for the abundance of daughter product s are the production of atomic 
carbon in the metastable *D state from CO ( Tozzi et al. 1998h and the high concentration of 
C02 + close to the nucleus seen in early IUE spectra (IFestou et al.lll982l ). The FUSE spectra 
provide additional evidence in the strong O I *D - 1 D° line at 1152.16 A, seen in Fig. HJ 
which also cannot be accounted for by solar resonance scattering alone. The lower state of 
this transition, the metastable *D state, is an abundant product of photodissociation of 
H2O, OH, and CO, and is the upper level of the oxygen "red" lines at 6300 a nd 6364 A, seen 



in ma ny comets and often used as a surrogate for the water production rate (IFeldman et al. 
20041 ). 



DISCUSSION 



4.1. Modeling 

The column density of the non-ther mal com p onent of CO was modeled assuming an 



H 2 CO p arent using the vecto ri al model of 



given by iHuebner et al.l (119921 ) . iMeier et al. 



Festoul (119811 ) with solar photodestruction rates 
Jl993h have also evaluated the photodissociation 



lifetimes (H2CO has two other dissociation channels, H + HCO, and H + H + CO) and find 
the rate into H? + CO to be 25% lower than th at of IHuebner et al.l while the total rate 
is 7% lower. iBockelee-Morvan & Crovisierl (119921) find a total photodissociation lifetime in 
close agreement with that of IHuebner et al.l The model assumes a nucleus source for the 
H 2 CO. Even if a distributed source were present (see Section |4~41 below) . the size of the 
FUSE aperture translates to a few thousand km projected on the sky and such a source can 
be ignored for our purposes. 

To determine production rates relative to water, we use values of Qh 2 o taken from the 
recent litera ture; values d erived from SOHO/SWAN measurements of cometary Lyman-o; 



emission by ICombi et al.l (12008 ), and those derived from ODIN submillimeter observations 
of H 2 directly by iBiver et al.l (120071 ). Where they o verlap, these two sets of values are 
in agreement to ~20-25%. We adopt the values from ICombi et al.l because of their finer 
sampling grid, allowing us to use data from the sam e day as the FUSE observations. We 
use the deconvolved model results from ICombi et al.l rather than the daily averages given. 



For C/2001 Q4 we use a value derived fro m near-simu ltaneous HST observations of OH 



emission, which is close to that reported by IBiver et al. 



For C/2001 A2 we note a significant difference from the value of Qn 2 o cited in the initial 
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reports of the FUSE observations by iFeldman et all (120021 ) and IWeaver et al.l (120021 ) . The 
earlier estimate was based on a visual outburst two days prior to the observation and a mod- 
eled production rate adjusted to mat ch the derived H? column density . The latter was based 
on fluorescence efficiencies derived by lKrasnopolsky fc Feldmanl (120011 ). With the availabilit y 
of detailed solar Lyman-/? fluxes and line shapes from SOHO/SUMER (ILemaire et al.ll2002l ). 
the heliocentric velocity dependent g-factors were recalculated and the H 2 column density 
revised downward by a fa ctor of about three, which is consistent with the water production 



rate for the date given by ICombi et al. 



The results are given in Table HI For the comets in which the non-thermal CO is 
detected, the production rates of H2CO relative to that for H2O are ge n erally in the range 
of values found from infrared and millimeter observations. iBiver et al.l (120061 ) reported on 
four comets observed at millimeter wavelengths including the three comets observed by 
FUSE at comparable epochs in 2001. For co mets C/2001 A2 and C/2000 WM1, our relative 
production rates are comparable to those of IBiver et al.l based on their model that assumes 
a distributed source with an effective parent lifetime about a factor of 2 l onger than the 
true photochemical lifetime. With an assumed H2CO parent, IBiver et al.l derive relative 
production rates 2 — 4 times smaller. 



Millimeter observations of comet C/2001 Q4 were made in mid-May 2004 by lMilam et al. 



(120061 ). Their results for H2CO, again assuming a distributed grain source, are in close 
agreement with the FUSE value. While infrared bands of H 2 CO near 3.6 [im have been 
detected in a number of recent com ets, the only dat a available f or comparison with the 
FUSE observations is of C/2001 A2 (jGibb et al.ll2007h . bibb et all observing 2 and 3 days 
prior to the FUSE observation, found a low value of 0.24% for the ratio of the production 
rate of H 2 CO relative to H 2 on July 9.5 followed by a decrease of a factor of 4 one day 
later. In contrast, we derive a much higher value for this ratio on July 12.5. The apparent 
variability of H 2 CO is discussed below in Section |4"72"1 

Modeling the H2 brightness is not as straightforward. The H2 molecules initially carry 
away most of the excess energy of dissociat ion as translationa l energy (~11 km s -1 ), but they 
are rapidly thermalized in the inner coma (ICombi et al.ll2004j ). The excited state populations 
will also vary with incident wavelength and be altered by collisions, so the actual number 
of H2 molecules in the v = 1, J = 3 level produced by solar photodissociation will have a 
large uncertainty. Nevertheless, as a consistency check we can evaluate the number of such 
molecules in the field-of-view needed to produce the observed 0.13 rayleigh brightness of the 
(1,1) Q(3) line at 1031.86 A observed in comet C/2001 Q4. The fluorescence efficiency, or 
g- factor, at the p eak of the sola r O VI line is calculated using the solar minimum irradiance 
recommended by I Warren! (120051 ) and the oscillator strengths and branching ratios taken from 
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the compilation of lAbgrall et al.l (ll993bl ). The derived g-factor is 3.8 x 10 7 s 1 molecule 



at 1 AU, so that a mean column density of 3.6 x 10 cm is sufficient to produce the 
observed emission rate. From Table (2J the column density of non-thermal CO in this comet 
i s 4.1 x 10 12 cm" 2 , a nd 78% of those are accompanied by H 2 rather than by two H atoms 
(IHuebner et al.lll992l ). Thus 11% of the total H 2 needs to be in v = 1, J = 3, which i s 



plausible based on the laboratory data of IZhang et al.l (120051 ) and IChambreau et al.l (120061 ). 
Given all of the uncertainties, this is good agreement and supports our proposal that H 2 CO is 
the source of both the non-thermal CO and the vibrationally excited H 2 detected by FUSE. 

The H 2 vibrational population may also be examined with the identification of solar 
Lyman-o; pumping of lines in the Lyman (1,2) band as described in Section [3731 The solar 
Lyman-a spectral flux is roughly two orders of magnitude larger than that of O VI A 103 1.9 
resulting in comparabl y larger fluore s cence efficiencies for these lines. Using solar minimum 
Lyman-q fluxes from iLemaire et al.l (120021 ) and molecular data for the Lyman bands from 
Abgrall et al.l (Il993al ). the g-factors for the (1,1) P(5) and R(6) lines can be calculated. The 
relative population of the v = 2, J = 5, needed to account for the observed 1161.8 A emission 
of 0.24 rayleighs is 2%. Thi s is consistent wit h the relative population of this level seen in 
the laboratory experiments (IZhang et al.ll2005l ). when we consider that these experiments are 
done at discrete wavelengths while in the cometary context photodissociation is produced 
by the wavelength integrated solar flux. 

We conclude that photodissociation of H 2 CO is sufficient to produce vibrationally ex- 
cited H 2 that accounts for the presence of Lyman-a excited fluorescence of H 2 in the coma. 
While the production rate of H 2 CO is ~0.5% relative to H 2 0, the solar photod issociation rate 
into H 2 is ~100 times higher for H 2 CO, depending on the solar activity level (IHuebner et al. 



19921 ). Thus both molecules contribute comparable amounts of H 2 to the coma. The branch- 
i ng of water photodiss ociation into H 2 + 0( X D) has also been studied in the laboratory 
( iSlanger fc Blackl 119821 ). but there are no available measurements of the resultant H 2 level 
populations. Although there has been extensive theoretical work on the H + OH channel of 
H 2 photodissociation (jSchinkelll993l ). there does not seem any devoted to the H 2 + 0( X D) 
channel, which is clearly warranted. O n the basis of the brightness of the Lyman-/3 pumped 
fluorescent lines, iFeldman et al.l (120021 ) suggested that solar photodissociation of water led 
to a thermal population of H 2 molecules in the v = 0, J = 1 level. 



4.2. Caveat on Temporal Variability 



We noted above the utility of the SOHO/SWAN data of ICombi et al.l (120081 ) to obtain 
daily values of the water production rate. For C/2001 A2 near the time of the FUSE ob- 
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servations, the comet exhibited several visual outbursts, particularly one on 2001 July 11.5, 
one d ay before the FUSE o bservations, which saw a 1.5 magnitude increase in visual bright- 
ness (ISekanina et al.ll2002l ). These out bursts were re flected in sudden increases in Qh 2 o of 
factors of 2-3 in one day in the data of ICombi et all The FUSE observations of this comet 
spanned 7.5 hours during which the CO emission decreased by a factor of two (Fig. [6]), likely 
a manifestation of the decrease in coma abundance following the outburst. Note that for 
this comet the FUSE 30" x 30" aperture projects to 6500 x 6500 km on the sky so that a 
molecule moving radially outward at 1 km s _1 would exit the field-of-view in less than one 
hour. Hence the observed variation represents variability of the CO source on this time scale. 
On the other hand, the column densities given in Table [2] are taken from the summed spectra 
(to improve the signal/noise ratio) and thus represent average values. Caution is therefore 
urged in deriving relative production rates from diverse observations, even if obtained on the 
same date, and in the comparison of results from different observers. 

Th e July 11.5 outbu rst makes it difficult to compare the derived H 2 CO production 
rates of iGibb et al.l (120071 ) on July 9.5 and 10.5, noted above, with the present results. From 
the limited data available, it is not possible to determine whether H 2 CO production was 
enhanced by the outburst or whether the increase inferred from the FUSE spe ctra is just 
the result of intrinsic variability responsible for the change seen by IGibb et al.l a few days 
earlier. The same holds for a comparison of ou r results with the CO pro duction rate derived 
from the same infrared spectrum on July 10.5 (IMagee-Sauer et al.ll2008l ) even though in this 
case the values of Qco are comparable. 

A similar caveat applies to comet C/2001 Q4 which FUSE observed over a 27 h period. 
The variation of the CO emission is shown in Fig. [7] and can be fit with a sinusoidal curve 
of period 17.0 h and 22% amplitude. Yet the CO production rate, derived from HST /STIS 
observations of the CO Fourth Positive system (ILupu et al.l 120071 ) . made near a minimum 
in the FUSE lightcurve, is consistent with the mean value derived from the FUSE spectra. 
In contrast, comet C/2001 WM1 showed little variation in CO emission over a 60 h period. 



For this comet, simultaneous HST/STIS observations were also obtained and ILupu et al. 
derived comparable CO production rates from both sets of spectra. 



4.3. Relative Abundance of CO2 

As noted in Section I3~l4"t for molecules with long lifetimes against solar photodissociation 
such as CO2, electron impact dissociation becomes an important destruction channel. Using 
the electron impact rates given in Table [31 however, would still require an abnormally large 
(> 20%) abundance of C0 2 to account for the derived "hot" CO column density. Since 
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there are no other suitable parents for CO other than CO2 and H2CO, it is likely that the 
model used to calculate the values given i n Table [3] signifi c antly underestimates the electron 
impact contribution to CO2 dissociation. I Cosby Helml (119921 ) note that the threshold for 
this process is at 7.4 eV whereas the ch annel leading to th e a 3 II state of CO requires 11.4 eV. 
At low electron energi es, the model of iKorosmezey et al.l (119871 ) does not match the data of 



Gringauz et al.l (119861 ). which appears to follow an E 1 distribution, so that there may be 



an additional large contribution to the dissociation rate from low energy electrons. There 
are also uncertainties in the radial and angular distribution of the electrons that make our 
estimate quantitatively unreliable. A more complete analysis of cometary photoelectrons, 
which is beyond the scope of the prese nt paper, wou l d enta il an analysis of CO Cameron band 
emission similar to that presented by I Weaver et al.l (11994J ) . Fortunately, relevant HST /STIS 
spectra do exist for comet C/2001 Q4 and will be treated in a separate publication. 



4.4. Distributed Source of CO 



The first evidence for a distributed source of CO in the coma of a comet came from the 
in situ Neutral Mass Spectrometer expe ri ment on the Giotto mission to comet lP/Halley 
(jEberhardt et al.l 119871 ; lEberhardtl Il999l ). iMeier et al.l (119931 ) also showed, from the same 
data set, that the H 2 CO also appeared to com e from an extend ed source and that the H 2 CO 
was sufficient to account for the extended CO (jEberhardtlll999l ) . Polyoxymethylene (POM), 
evaporated fro m organic materia l on refractory grains, has been postulated as the source 
of this H 2 CO (jCottin et al.l 12004 ). In recent years, as high resolution spectroscopic instru- 
mentation in the near infrared has made the observation of CO in comets almost routine, 
considerable effort has been expended in trying to determine if such distributed sources are 
common or unique to certain comets (IBockelee-Morvan et al.l 120041 ) . Identification of a dis- 
tributed emission source is based on long-slit spectroscopy in which the spatia l profile of the 



emiss ion deviates from that expected for radial outflow (IDiSanti et al.l l200ll ; iBrooke et al. 
20031 ). For CO this can be done in both the infrared and the ultraviolet. 



This approach is applicable only to the brightest comets because the spatial pixels along 
the slit must be small enough to provide adequate spatial resolution. For such comets, such 
as C/1995 01 (Hale-Bopp), the CO column densities along the line-of-sight (to both the 
Sun and to the observer) are sufficiently large that optical depth effects tend to produce a 
flattening of the spatial profile near the nucl eus that mimics th e effect of a distributed source 
(th e ultraviolet c ase has been discussed by iLupu et al.l 120071 ) . This concern was addressed 
by IDiSanti et al.l in their analysis of CO in comet Hale-Bopp. They found a very high 
production rate of CO relative to that of water, 24.1%, with half being from the nucleus and 
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the other half from a distributed source. In the same comet, iBrooke et all (120031 ) found an 
even higher relative abundanc e of CO, 37-41%, with 90% coming from a distributed source. 
Bockelee-Morvan et all (120051 ) questioned both of these results su gge sting that the y were the 
effect of a large infrared opacity, even though both lDiSanti et all and lBrooke et all accounted 
for opacity effects. In any case, if such a distributed source were H 2 CO it would imply a n 
abundance in this comet an order of magnitude higher than that cited by lBiver et all (120061 ). 



5. CONCLUSION 

Spectroscopic observations of four comets by FUSE have opened a new window for the 
investigation of the physics and chemistry of the cometary coma. The high spectral resolution 
of FUSE has made it possible to separate the multiple components of the CO emission, while 
the high sensitivity has allowed the detection of fluorescent H 2 emission pumped by solar 
O VI radiation. Both of these can be understood in terms of the dissociation processes of 
H 2 CO and C0 2 . Despite the small sample, it is clear that previous ultraviolet spectroscopic 
observations of CO in comets including those from sounding rockets, IUE and HST (e.g., 



Feldman et all 119971 ) probably included these additional components of CO that were not 
resolved, either spectrally or spatially, from the native ( "cold" ) CO with the instrumentation 
used. Photoelectron excitation contributes ~10% of the observed emission of the CO C — X 
(0,0) band and needs to be accounted for in deriving CO production rates from FUSE 
observations. 
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Table 1: Comet observation parameters and observed CO brightnesses. 



Comet 


C/1999 Tl 


C/2001 A2 


C/2000 WM1 


C/2001 Q4 




(McNaught-Hartley) 


(LINEAR) 


(LINEAR) 


(NEAT) 


Date 


2001 Feb 3 


2001 Jul 12 


2001 Dec 7-10 


2004 Apr 24 


Start Time (UT) 


17:43 


13:48 


08:48 


00:40 


Exposure Duration (s) 


15,614 


15,063 


36,467 


68,282 


r (AU) 


1.43 


1.20 


1.12 


1.03 


A (AU) 


1.29 


0.30 


0.34 


0.51 


f (km s _1 ) 


14.9 


22.8 


-28.3 


-10.8 


F10.7 (sfu) 


160 


140 


220 


115 


C -X (0,0) (rayleighs) 


1.46 ±0.03 


2.47 ±0.04 


0.74 ±0.02 


6.29 ±0.03 


. . .P(40) peak (rayleighs) 


< 0.02 


0.16 ±0.03 


0.10 ±0.01 


0.16 ±0.02 


B - X (0,0) (rayleighs) 


0.19 ±0.02 


0.71 ±0.03 


0.16 ±0.01 


1.51 ±0.02 



Table 2: Derived CO column densities and rotational temperatures in the FUSE 30" x 30" 
aperture. Errors are from the la contours in a three parameter fit. 



Comet 




C/1999 Tl 


C/2001 A2 


C/2000 WM1 


C/2001 Q4 






(McNaught-Hartley) 


(LINEAR) 


(LINEAR) 


(NEAT) 


A^ co i d (cm 


- 2 ) 


(3.5 ±0.2) x 10 13 


(2.7 ±0.15) x 10 13 


(6.0 ±0.4) x 10 12 


(1.55 ±0.05) x 10 14 


A hot (cm" 


- 2 ) 


< 1.3 x 10 12 


(1.10 ±0.06) x 10 13 


(2.5 ±0.5) x 10 12 


(2.0 ±0.1) x 10 13 


A^non— thermal (cm ) 


< 0.8 x 10 12 


(4.5 ±0.7) x 10 12 


(2.5 ±0.3) x 10 12 


(3.3 ±0.3) x 10 12 


T cold (K) 




59 ±5 


56 ±7 


75 ±10 


58 ±2 


T hot (K) 




500 


500 


500 


600 



a Assumes an equal contribution from the R(40) peak not observed. Uncertainties are la in the counting 
statistics. 
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Table 3: Derived effective electron impact excitation rates (</), in molecule 1 s 1 , averaged 
over the FUSE aperture. 



Comet 


C/f999 Tf 


C/200f A2 


C/2000 WM1 


C/2001 Q4 




(McNaught-Hartley) 


(LINEAR) 


(LINEAR) 


(NEAT) 


B-X (0,0) a 


9.0 x 10~ 9 


2.4 x I(T 8 


2.1 x 1(T 8 


7.4 x 1(T 9 


C-X (0,0) a 


5.8 x lO- 9 


1.6 x iir 8 


i.3 x iir 8 


4.7 x 10~ 9 


a — X (all bands) b 


2.6 x f(T 7 


7.1 x 10~ 7 


6.1 x 10~ 7 


2.2 x 10~ 7 


total dissociation 15 


f.f x f(T 6 


2.9 x iir 6 


2.4 x 1(T 6 


8.6 x 1(T 7 



"electron impact on CO at 1 AU. For comparison, the optically thin fluorescence efficiencies for the B — X (0,0) 
and C — X (0,0) bands at 1 AU are 2.3 x 10~ 9 and 1.0 x 10~ 7 molecule -1 s _1 , respectively, for solar maximum 
conditions. 

b electron impact on CO2 at 1 AU. 



Table 4: Modeled production rates and production rate ratios. 



Comet 


C/1999 Tl 


C/2001 A2 


C/2000 WM1 


C/2001 Q4 




(McNaught-Hartley) 


(LINEAR) 


(LINEAR) 


(NEAT) 


Qco (s" 1 )" 


7.8 x 10 27 


1.33 x 10 27 


3.7 x 10 26 


1.45 x 10 28 


<5h 2 co (s _1 ) 


< 3.7 x 10 26 


7.6 x 10 26 


4.6 x 10 26 


8.1 x 10 26 


Qu 2 o (s- 1 ) 


7.6 x 10 28 


6.0 x 10 28 


8.8 x 10 28 


2.0 x 10 29 


Qco/Qh 2 o (%) a 


10.3 


2.2 


0.41 


7.2 


QH2C0/QH20 (%) 


<0.5 


1.3 


0.52 


0.40 



"corrected for photoelectron impact contribution. 
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FIGURE CAPTIONS 



Fig. 1.— CO C-X (0,0) band in comets McNaught-Hartley, C/2001 A2, C/2000 WM1 and 
C/2001 Q4. The model fit is shown together with the "cold" (dotted) and "hot" (dashed) 
components. P(40) indicates the position of the P-branch transitions originating from J 
levels near the peak of the hot, non-thermal ro tational distribution seen in laboratory pho- 



todissociation experiments (Ivan Zee et al. 



19931 ) . The corresponding R-branch lines lie below 



the short wavelength cut-off of the LiF2a detector. For C/2001 Q4, the gray line represents 
the vertical scale expanded by a factor of 10 to enhance the visibility of the non-thermal 
peak. 

Fig. 2. — O VI pumped H2 lines in the FUSE spectrum of comet C/2001 Q4. The wavelength 
scale is in the comet's rest frame. The exposure time for the LiFla channel (left panel) was 
67,360 s, while the other two panels include the sum of LiFlb and LiF2a with a total 
exposure time of 137,530 s. The feature indicated by (*) is identified as solar Lyman-a 
pumped fluorescence of H2 in the Lyman (1,1)P(5) and R(6) lines. The Lyman (6,3)P(1) 
line is pumped by solar Lyman-/?. 

Fig. 3. — Solar O VI line profile derived from an instrumentally scattered FUSE spectrum. 
The position of the Werner (1,1) Q(3) line is shown as the dashed line. The heliocentric 
Doppler shifts of the four comets studied are indicated as (a) C/2001 A2, (b) McNaught- 
Hartley, (c) C/2001 Q4, (d) C/2000 WM1. 

Fig. 4. — Same as Fig. [1] for the CO B — X (0,0) band. The predicted solar fluorescence 
brightness is shown for both the "cold" (dotted) and "hot" (dashed) components. The shape 
of the neighboring O I A1152.15 feature reflects the spatial distribution of the emitter in the 
FUSE aperture. 



Fig. 5. — Cross sections for electron impact excitation (from lShirai et al.ll200ll ). The curves 
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labeled (B) and (C) are for excitation of CO to produce B — X (0,0) and C — X (0,0) band 
photons, respectively, while (a) represents electron impact dissociation of CO2 into the a 3 I1 
state of CO. The dashed cur ve represents a theoretical photoelectron energy distribution 



(from iKorosmezey et al.l 119871 ) . 



Fig. 6. — Temporal variability of the CO C — X and B — X bands in comet C/2001 A2. 
The horizontal bars indicate the individual exposures. The visual brightness deduced from 
the FUSE Fine Error Sensor is also shown. 

Fig. 7. — Temporal variability of the CO C — X band in comet C/2001 Q4. The origin of 
time is UT 2004 April 24 00:40:32. The solid line is a sinusoidal fit with a period of 17.0 h 
and an amplitude of 22%. 
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Fig. 1.— CO C-X (0,0) band in comets McNaught-Hartley, C/2001 A2, C/2000 WM1 and 
C/2001 Q4. The model fit is shown together with the "cold" (dotted) and "hot" (dashed) 
components. P(40) indicates the position of the P-branch transitions originating from J 
levels near the peak of the hot, non-thermal ro tational distribution seen in laboratory pho- 
todissociation experiments (Ivan Zee et al.lll993l ). The corresponding R-branch lines lie below 
the short wavelength cut-off of the LiF2a detector. For C/2001 Q4, the gray line represents 
the vertical scale expanded by a factor of 10 to enhance the visibility of the non-thermal 
peak. 
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Fig. 2. — O VI pumped H 2 lines in the FUSE spectrum of comet C/2001 Q4. The wavelength 
scale is in the comet's rest frame. The exposure time for the LiFla channel (left panel) was 
67,360 s, while the other two panels include the sum of LiFlb and LiF2a with a total 
exposure time of 137,530 s. The feature indicated by (*) is identified as solar Lyman-a 
pumped fluorescence of H 2 in the Lyman (1,1)P(5) and R(6) lines. The Lyman (6,3)P(1) 
line is pumped by solar Lyman-/?. 
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Fig. 3. — Solar O VI line profile derived from an instrumentally scattered FUSE spectrum. 
The position of the Werner (1,1) Q(3) line is shown as the dashed line. The heliocentric 
Doppler shifts of the four comets studied are indicated as (a) C/2001 A2, (b) McNaught- 
Hartley, (c) C/2001 Q4, (d) C/2000 WM1. 
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Fig. 4. — Same as Fig. [T] for the CO B — X (0,0) band. The predicted solar fluorescence 
brightness is shown for both the "cold" (dotted) and "hot" (dashed) components. The shape 
of the neighboring O I A1152.15 feature reflects the spatial distribution of the emitter in the 
FUSE aperture. 
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Fig. 5. — Cross sections for electron impact excitation (from 
labeled (B) and (C) are for excitation of CO to produce B — X (0,0) and C — X (0,0) band 
photons, respectively, while (a) represents electron impact dissociation of CO2 into the a 3 fl 
state of CO. The dashed cur ve represents a theoretical photoelectron energy distribution 
(from iKorosmezey et al.l 119871 ) . 
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Fig. 6.— Temporal variability of the CO C - X and B-X bands in comet C/2001 A2. 
The horizontal bars indicate the individual exposures. The visual brightness deduced from 
the FUSE Fine Error Sensor is also shown. 
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Fig. 7. — Temporal variability of the CO C — X band in comet C/2001 Q4. The origin of 
time is UT 2004 April 24 00:40:32. The solid line is a sinusoidal fit with a period of 17.0 h 
and an amplitude of 22%. 



